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The  pressure  gradients  in  the  electrodes  of  a  solid  oxide  fuel  cell  (SOFC)  are  frequently  neglected  without 
any  justification  in  calculating  the  concentration  overpotentials  of  the  SOFC  electrodes  in  modeling  stud¬ 
ies.  In  this  short  communication,  a  comparative  study  has  been  conducted  to  study  the  effect  of  pressure 
gradients  on  mass  transfer  and  the  resulting  concentration  overpotentials  of  an  SOFC  running  on  methane 
(CH4)  fuel.  It  is  found  that  the  pressure  gradients  in  both  anode  and  cathode  are  significant  in  the  fuel  cell 
electrochemical  activities.  Neglecting  the  anode  pressure  gradient  in  the  calculation  can  lead  to  under¬ 
estimation  of  the  concentration  overpotential  by  about  20%  at  a  typical  current  density  of  5000  Am  2 
and  at  a  temperature  of  1073  K.  The  deviation  can  be  even  larger  at  a  higher  temperature.  At  the  cathode, 
neglecting  the  pressure  gradient  can  result  in  overestimation  of  the  concentration  overpotential  by  about 
10%  under  typical  working  conditions. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  promise  efficient  electricity  gen¬ 
eration  for  stationary  applications  at  an  efficiency  greater  than 
conventional  heat  engines  [1].  Operating  at  high  temperatures 
(1073-1273 1<),  SOFCs  have  a  number  of  advantages  over  low- 
temperature  fuel  cells  such  as:  (1)  electrochemical  reactions  are 
fast  at  high  temperatures,  leading  to  low  activation  overpotential 
and  effective  use  of  low-cost  catalyst,  such  as  Ni;  (2)  the  electrolyte 
at  a  high  temperature  has  a  high  ion  conductivity,  thus  the  ohmic 
overpotential  of  the  electrolyte  can  be  minimized;  (3)  high  operat¬ 
ing  temperature  enables  direct  internal  reforming  of  hydrocarbon 
fuels  in  SOFCs,  thus  SOFCs  can  utilize  a  variety  of  fuels,  including 
CO,  which  is  an  unwanted  poisonous  gas  for  low-temperature  fuel 
cells;  (4)  waste  heat  from  SOFCs  is  of  high  quality  and  can  be  recov¬ 
ered  by  using  a  bottoming  cycle  to  increase  the  system  efficiency 
[2]. 

However,  the  high  operating  temperature  limits  the  choice  of 
materials  used  for  SOFCs.  The  adverse  thermal  expansion  mismatch 
of  the  SOFC  components  and  catalyst  sintering  may  also  occur  at  a 
high  temperature.  In  order  to  resolve  these  material  and  long-term 
stability  problems,  great  research  efforts  have  been  done  to  reduce 
the  operating  temperature  of  SOFCs  to  about  773  K  [3].  The  reduc¬ 
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tion  in  operating  temperature  decreases  the  electrode  activity  and 
electrolyte  ionic  conductivity,  which  in  turn  leads  to  lower  SOFC 
performance.  Accordingly,  the  development  of  electrode  with  high 
electrochemical  activity  and  electrolyte  with  high  ionic  conduc¬ 
tivity  are  desired.  An  alternative  way  to  reduce  the  considerable 
ohmic  overpotential  is  to  fabricate  thin  film  electrolyte.  Presently, 
thin  electrolyte  with  a  thickness  of  only  10  pim  can  be  fabricated 
by  a  number  of  methods,  such  as  tape  casting  and  co-firing  [4], 
sol-gel  method  [5],  and  polymer  assisted  combustion  method  [6]. 
In  order  to  provide  a  strong  mechanical  support  for  the  SOFC,  either 
the  anode  or  the  cathode  must  be  thick  enough  (>500  pim).  Anode 
support  is  the  most  favorable  configuration  for  SOFCs  as  the  cell 
performance  is  generally  much  better  than  the  cathode-supported 
SOFC  [7-9]. 

In  an  anode-supported  SOFC,  accurate  calculation  of  the 
gas  composition  at  the  anode-electrolyte  interface  is  extremely 
important  as  it  can  strongly  affect  not  only  the  concentration 
overpotentials  but  also  the  activation  overpotentials,  due  to  the 
dependence  of  the  electrochemical  reaction  kinetics  on  the  con¬ 
centrations  of  reactant  and  product  species  [8,10].  In  particular, 
accurate  calculation  of  the  concentration  overpotential  is  very 
important  for  design  optimization  of  an  anode-supported  SOFC 
working  at  high  fuel  utilizations,  as  “a  small  error  in  the  concentra¬ 
tion  overpotential  calculation  may  cause  a  dramatic  change  to  its 
design  performance”  [11]. 

The  concentration  overpotentials  of  an  SOFC  can  be  modeled 
by  Fick  model  (FM),  Dusty  Gas  model  (DGM)  or  Stefan-Maxwell 
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model  (SMM).  Among  the  three  models,  DGM  can  predict  the  mass 
transfer  in  porous  SOFC  electrodes  more  accurately  as  it  takes 
into  account  Knudsen  diffusion  phenomena  and  eliminates  the 
assumption  of  equal-molar  counter  diffusion  in  FM,  which  becomes 
invalid  if  the  molecular  weights  of  the  diffusing  gas  species  differ 
considerably  [11  ].  In  addition  to  diffusion,  another  possible  mecha¬ 
nism  governing  the  multi-component  mass  transfer  is  permeation, 
which  is  driven  by  pressure  gradient  inside  the  porous  electrodes. 
This  pressure  gradient  can  be  produced  by  non-equal-molar  diffu¬ 
sion  of  gas  species  [12],  or  different  molar  consumption/production 
of  reactants/products.  For  example,  in  the  porous  cathode,  02  gas 
is  consumed  without  any  gas  product;  thus,  a  pressure  gradient 
will  be  established  across  the  porous  cathode  layer.  In  the  porous 
anode  of  a  H2  fed  SOFC,  the  pressure  gradient  occurs  due  to  differ¬ 
ent  diffusion  rates  of  H2  and  H20.  When  hydrocarbon  fuels,  such 
as  methane,  are  used  in  SOFCs,  internal  methane  steam  reform¬ 
ing  (MSR)  and  water  gas  shift  (WGS)  reactions  can  be  initiated  in 
the  porous  anode,  which  in  turn  can  lead  to  a  non-uniform  dis¬ 
tribution  of  pressure  in  the  SOFC  anode.  However,  most  studies 
in  literature  totally  neglected  the  effect  of  pressure  gradient  on 
concentration  overpotentials  of  SOFC  electrodes  without  any  justi¬ 
fication  even  for  SOFC  running  on  hydrocarbon  fuels  with  MSR  and 
WGS  [7,10,13-46].  Therefore,  it  is  still  not  clear  how  and  to  what 
extent  the  existence  of  pressure  gradient  affects  the  concentration 
overpotential  calculation.  This  short  communication  is  purposely 
designed  to  study  the  effect  of  pressure  gradients  on  SOFC  con¬ 
centration  overpotentials.  In  this  study,  methane  is  considered  as  a 
model  hydrocarbon  fuel  for  the  SOFC. 

2.  The  model 


2.1.  Concentration  overpotentials 


Concentration  overpotential  reflects  the  resistance  of  the 
porous  electrode  structure  to  the  transport  of  reacting  species 
to  and  products  from  the  triple  phase  boundary  (TPB)  at  the 
electrode-electrolyte  interface,  where  the  electrochemical  reac¬ 
tions  take  place.  For  SOFCs  using  pure  H2  fuel,  or  when  H2  is  the  only 
electrochemically  reacting  fuel,  the  concentration  overpotential  at 
the  anode  (^conc.a)  and  cathode  (^conc.c)  can  be  expressed  as 
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P  cone,  a  —  2p  n 
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sidered.  DGM  takes  into  account  the  molecular  diffusion,  Knudsen 
diffusion,  and  permeation  for  the  multi-component  mass  transfer 
in  the  porous  electrodes, 
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where  Afz-,yz-,  and  are  the  flux,  molar  fraction  and  effective  Knud¬ 
sen  diffusion  coefficient  of  species  i;  D?.ff  is  the  effective  binary 
diffusion  coefficient  of  species  i  and  j;  P  is  pressure;  x  is  the  depth 
measured  from  electrode  surface;  fi  is  the  viscosity  of  the  gas  mix¬ 
tures;  and  B0  is  the  permeability  of  the  porous  electrode,  which  can 
be  calculated  by  the  Kozeny-Carman  relationship  [47], 
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The  effective  diffusion  coefficients  (D?f£,  D9.ff)  can  be  determined 
by  the  following  equations  [48], 
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where  £ ,  £,  and  rp  are  the  porosity,  tortuosity  and  mean  pore  radius 
of  the  electrodes;  and  M*  is  the  molecular  weight  of  species  i;  T2D 
is  a  dimensionless  diffusion  collision  integral  and  critj  is  the  mean 
characteristic  length  of  species  i  and  j.  The  values  of  T2D  and  a can 
be  obtained  from  literature  [49]. 

The  pressure  gradient  (dP/dx)  can  be  evaluated  with  a  method 
developed  by  Zhu  and  Kee  [50]: 

d p  =  -EL  iW/pff) 

dx  (\/RT)  +  (B0P/RTti)Y,l1iyi/Dfk) 

Eq.  (17)  comes  from  summation  of  Eq.  (3)  over  all  species  i. 

If  the  effect  of  pressure  gradient  is  neglected,  then  Eq.  (3)  is 
simplified  to 
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which  is  commonly  used  in  literature. 


2.3.  Methane  steam  reforming  and  water  gas  shift  reactions 


where  R  is  the  universal  gas  constant;  T  is  the  working  tempera¬ 
ture;  Pis  the  Faraday  constant;  P2  is  the  partial  pressure  of  species 
i  (i  refers  to  H2,  H20,  and  02  for  respective  pressure);  the  super¬ 
scripts  S  and  I  indicate  the  partial  pressure  at  the  electrode  surface 
and  the  electrode-electrolyte  interface,  respectively.  The  surface 
partial  pressures  (P^2,  P^  0»  ^o2 )  are  giyen  as  known  input  param¬ 
eters  while  the  interface  partial  pressures  (Pf^ ,  P{j  0,  anc*  ^o2 ) nee^ 
to  be  determined  by  analyzing  the  mass  transfer  inside  the  porous 
electrodes. 

2.2.  Multi-component  mass  transfer  in  porous  electrodes 

The  mass  transfer  in  porous  SOFC  electrodes  is  governed  by  two 
mechanisms:  (1 )  diffusion  due  to  concentration  gradients  of  react¬ 
ing  species  and  (2)  permeation  driven  by  pressure  gradient.  As  the 
mean  pore  size  of  SOFC  electrode  is  comparable  to  the  mean  free 
path  of  the  reacting  gases,  the  Knudsen  diffusion  effect  must  be  con- 


In  practice,  an  SOFC  can  be  fed  with  H2  or  hydrocarbon  fuels. 
Here,  CH4  is  used  as  a  model  hydrocarbon  fuel  for  SOFC.  The  major 
chemical  reactions  involved  are  methane  steam  reforming  (MSR) 
and  water  gas  shift  reactions, 

CH4  +  H20  CO  +  3H2  (9) 

CO  +  H20  C02  +  H2  (10) 

It  should  be  mentioned  that  both  H2  and  CO  could  be  electro¬ 
chemically  oxidized  at  the  TPB.  However,  it  is  legitimate  to  neglect 
the  electrochemical  oxidation  of  CO  because:  (1)  the  rate  of  CO 
electrochemical  oxidation  is  much  lower  than  the  electrochemi¬ 
cal  oxidation  of  H2  and  (2)  the  rate  of  CO  electrochemical  oxidation 
is  much  lower  than  the  chemical  oxidation  of  CO  by  WGS  [51,52]. 

The  reaction  rates  (mol  m-3  s-1 )  of  MSR  (Rmsr)  and  WGS  (Rwgs) 
in  the  porous  anode  can  be,  respectively,  written  as  [53]. 

Rmsr  =  /<mSRPch4Ph2o  -  RmsrPC0^2  )3  (H ) 
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and 


flwGS  =  fcwGsPcoPH20  “  ^wgS^C02Ph2  O2) 

where  k  are  rate  constants  for  MSR  and  WGS  (molnrr3  Pa-2  s-1); 
the  superscripts  +  and  -  refer  to  the  forward  and  backward  reac¬ 
tions,  respectively. 

In  a  steady  state,  the  transport  of  each  participating  component 
is  determined  by  the  local  conservation  of  mass: 
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i  and  x  is  the  depth 


measured  from  electrode  surface. 

The  governing  equations  (Eqs.  (3),  (7),  (8)  and  (11 )— (17))  can  be 
solved  numerically  to  obtain  the  partial  pressure  of  gas/product 
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Fig.  1.  Comparison  between  the  present  simulation  results  and  literature  data  [25] 
for  model  validation. 


species  at  the  electrode-electrolyte  interface.  It  should  be  men¬ 
tioned  that  in  the  calculation  procedure,  Eq.  (3)  should  be  applied 
to  all  transporting  species  to  model  the  transport  of  all  gas  species. 
The  finite  difference  method  can  be  used  to  discretize  the  gov¬ 
erning  equations,  which  can  be  solved  with  an  iterative  scheme. 


Fig.  2.  Distributions  of  pressure,  gas  composition  and  MSR  reaction  rate  in  the  anode  of  an  SOFC  running  on  CH4  with  internal  MSR  and  WGS:  (a)  pressure;  (b)  molar  fraction 
of  H2 ;  (c)  molar  fraction  of  H20;  and  (d)  reaction  rate  of  MSR. 
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Table  1 

Input  parameters  used  in  the  mathematical  model 


Parameter 

Value 

Temperature,  T(K) 

1073 

Pressure  at  the  electrode  surface,  P  (atm) 

1.0 

Fuel  composition  (molar  fraction)  at  the  anode  inlet  [25,52] 

ch4 

17.1% 

CO 

2.9% 

H20 

49.3% 

h2 

26.3% 

co2 

4.4% 

Electrode  porosity 

0.4 

Electrode  tortuosity 

5.0 

Electrode  pore  radius  (p,m) 

0.5 

Anode  thickness,  da  (p,m) 

500 

Cathode  thickness,  dc  (|Jim) 

50 
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After  obtaining  the  interfacial  partial  pressures  of  the  reacting  gas 
species,  the  concentration  overpotentials  of  the  anode  and  cathode 
can  thus  be  calculated  by  Eqs.  (1 )  and  (2),  respectively. 

3.  Results  and  discussions 

The  electrochemical  model  presented  in  Section  2  was  validated 
in  the  previous  studies  [54,55].  A  comparison  between  the  present 
simulation  results  and  data  from  literature  is  illustrated  in  Fig.  1. 
It  can  be  seen  that  the  present  results  agree  well  with  Ivanov’s 
[25]  data,  thus  the  model  presented  in  this  paper  is  reliable  for 
more  analyses.  In  this  section,  the  effect  of  pressure  gradient  on 
the  concentration  overpotentials  is  presented  and  discussed.  In  the 
following  analyses,  two  cases  are  considered:  (1)  with  pressure 
gradient  inside  the  porous  electrodes,  and  (2)  without  pressure  gra¬ 
dient  inside  the  porous  electrodes.  The  typical  values  of  the  model 
parameters  are  summarized  in  Table  1. 

3.1.  Anode  of  an  SOFC  running  on  CH4 

Fig.  2  shows  the  pressure,  gas  composition  and  MSR  reaction 
rate  in  the  anode  of  an  SOFC  running  on  CH4  with  internal  MSR  and 
WGS.  The  pressure  in  the  anode  is  found  to  increase  significantly 
along  the  anode  depth  (Fig.  2(a)).  This  is  because  (1 )  the  total  molar 
number  of  the  MSR  (Eq.  (9))  product  is  twice  as  many  as  that  of  the 
reactants  and  (2)  the  diffusion  of  H2  is  faster  than  that  of  H20  due 
to  its  smaller  molecular  weight  (non-equal-molar  diffusion)  [12]. 
As  this  pressure  gradient  adversely  impedes  the  transport  of  H2, 
the  H2  molar  fraction  is  found  lower  than  the  case  of  neglecting  the 
pressure  gradient  (Fig.  2(b)).  On  the  other  hand,  as  the  pressure  gra¬ 
dient  favors  the  transport  of  H20,  the  molar  fraction  of  H20  is  larger 
if  the  pressure  gradient  in  the  anode  is  considered  (Fig.  2(c)).  Due  to 
smaller  H2  molar  fraction  and  larger  H20  molar  fraction,  higher  rate 
of  MSR  is  observed  with  pressure  gradient  considered  (Fig.  2(d)). 
The  rate  of  MSR  initially  decreases  with  increasing  anode  depth  due 
to  the  decreased  molar  fraction  of  CH4  while  it  tends  to  increase  in 
a  layer  near  the  anode-electrolyte  interface  where  the  molar  frac¬ 
tion  of  H20  is  high.  This  explains  why  the  pressure  increases  at  a 
slightly  smaller  rate  in  deeper  layer  of  the  anode  (may  not  be  easily 
observed  in  Fig.  2(a)).  With  an  increase  in  current  density,  the  molar 
fraction  of  H2  decreases  while  the  molar  fraction  of  H20  increases 
due  to  the  enhanced  electrochemical  oxidation  reaction  (Figs.  2(b) 
and  (c)).  The  reduction  in  H2  molar  fraction  and  increase  in  H20 
molar  fraction  can  increase  the  rate  of  MSR  reaction  (Fig.  2(d)), 
which  in  turn  increases  the  pressure  gradient  in  the  porous  anode 
(Fig.  2(a)). 

When  the  pressure  gradient  inside  the  anode  is  considered,  the 
anode  concentration  overpotential  is  found  higher  than  the  case  if 


Fig.  3.  Effect  of  pressure  gradient  on  predicting  the  anode  concentration  overpo¬ 
tentials  of  an  SOFC  running  on  CH4:  (a)  varying  current  density:  and  (b)  varying 
temperature. 

the  pressure  inside  the  anode  is  assumed  to  be  invariant  (Fig.  3(a)). 
The  concentration  overpotential  is  actually  the  combined  effects 
of  mass  transfer  resistance  and  direct  internal  MSR  and  WGS.  At 
a  low  current  density,  the  molar  fraction  of  H2  at  the  TPB  can  be 
higher  than  that  at  the  anode  surface  due  to  H2  production  through 
MSR  and  WGS.  Therefore,  the  anode  concentration  overpotential 
is  found  negative  at  low  current  density  (Fig.  3(a)).  More  impor¬ 
tantly,  the  difference  in  concentration  overpotential  between  the 
two  cases  is  found  significant.  For  example,  at  a  typical  current 
density  of  5000  Am-2,  the  anode  concentration  overpotential  con¬ 
sidering  the  pressure  gradient  is  0.053  V,  while  the  value  is  0.036  V 
if  the  pressure  gradient  in  the  anode  is  totally  neglected.  In  addition 
to  current  density,  fuel  utilization  is  another  useful  parameter  for 
characterizing  the  SOFC  performance.  For  an  SOFC  working  under 
given  conditions  (temperature,  pressure,  flow  rate,  etc.),  the  fuel 
utilization  is  proportional  to  the  current  density.  Therefore,  the 
above  analyses  on  current  density  also  reveal  the  fact  that  high  fuel 
utilization  can  cause  high  concentration  overpotential. 

Fig.  3(b)  shows  the  effect  of  temperature  on  the  anode  concen¬ 
tration  overpotentials  at  a  typical  current  density  of  5000  Am-2. 
For  both  cases,  the  anode  concentration  overpotentials  are  found 
to  decrease  with  increasing  temperature.  This  is  different  from 
our  previous  analysis  on  H2  fed  SOFC,  in  which  the  concentration 
overpotential  is  found  to  increase  with  increasing  temperature  [8]. 
When  CH4  is  used  as  a  fuel  for  SOFC,  higher  temperature  favors 
MSR  and  WGS,  which  contribute  to  more  H2  production  inside  the 
porous  anode.  This  in  turn  increases  the  H2  molar  fraction  and  thus 
leads  to  lower  anode  concentration  overpotential.  In  addition,  the 
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difference  in  concentration  overpotential  between  the  two  cases 
becomes  more  significant  at  higher  temperature  (Fig.  3(b)).  Since 
MSR  is  favored  at  elevated  temperature,  the  pressure  gradient  is 
higher,  which  in  turn  leads  to  more  significant  effect  on  anode 
concentration  overpotential. 

It  should  be  mentioned  that  the  present  study  does  not  include 
energy  equation  and  assumes  uniform  temperature  distribution 
along  the  depth  of  porous  anode.  Actually,  the  endothermic  nature 
of  the  MSR  can  lower  the  temperature  of  the  porous  anode  and 
the  transporting  gas  species,  which  in  turn  can  reduce  the  rate  of 
MSR  for  H2  production.  As  a  result,  the  pressure  gradient  can  be 
less  significant  and  the  molar  fraction  of  H2  can  be  lower.  How¬ 
ever,  based  on  previous  analysis  on  temperature  distribution  in  the 
porous  anode  of  SOFC  running  on  methane,  the  temperature  vari¬ 
ation  in  the  anode  depth  is  much  less  than  10  °C  [52,56].  Thus,  the 
endothermic  nature  of  MSR  only  has  very  small  effect  on  distri¬ 
butions  of  pressure  gradient  and  gas  composition  along  the  anode 
depth.  Nevertheless,  the  endothermic  MSR  can  cause  considerable 
temperature  variation  along  the  fuel  flow  channel  because  of  large 
variation  in  gas  composition  along  the  channel  [52].  Therefore,  in 
2D  or  3D  models,  heat  transfer  model  should  be  included  to  account 
for  the  effect  of  reaction  heat. 

From  the  above  analyses,  it  can  be  seen  that  the  difference 
between  the  two  cases  is  significant  in  the  typical  range  of  cur¬ 
rent  density,  especially  at  a  high  temperature.  Neglecting  the  effect 
of  pressure  gradient  can  underestimate  the  concentration  overpo¬ 


tential  by  about  20%  at  a  typical  current  density  of  5000  A  m-2  (the 
anode  concentration  overpotential  starts  from  a  negative  value).  It 
is  thus  recommended  that  the  pressure  gradient  should  be  included 
in  modeling  the  anode  concentration  overpotential  of  an  SOFC  run¬ 
ning  on  hydrocarbon  fuels. 

3.2.  Cathode  of  an  SOFC 

In  this  section,  the  mass  transfer  and  concentration  overpoten¬ 
tial  of  the  cathode  are  investigated.  Different  from  the  anode,  no 
chemical  reactions  are  involved  in  the  cathode.  Taking  air  as  the 
oxidant,  the  pressure  gradient  in  the  cathode  can  be  expressed  as, 

dPc  =  -((N02/<,l<)-(^N2/<k)) 

dx  ((\/RT)  +  (B0Pc/RT^)mo2/Df2M)  +  {yn2/Df2k)) 

It  is  noted  that  Nq2  is  equal  to  J/4F  throughout  the  cathode  layer 
while  NN2  is  equal  to  zero.  Considering  that  there  is  small  difference 
in  molecular  weights  between  of  02  and  N2,  k  can  be  approx¬ 
imated  by  k.  Assuming  constant  viscosity  (/x)  in  the  porous 
cathode,  Eq.  (18)  can  be  re-organized  to 

g.  T  „9, 

dx  (1/RT)  +  (BoPc/RT[iDq2  k) 

In  Eq.  (19),  only  Pc  is  a  function  of  cathode  depth  x,  while  all  other 
parameters  are  independent  of  x.  With  the  boundary  condition 
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of  Pc|x_0  =PS  =  101310.0  Pa,  the  analytical  solution  of  the  above 
pressure  equation  can  be  derived  as 


do2mRT/1 

(  1  .  BoPcS  N 

2 

\  2BoNo2x  1 

Bo 

\ 

VRT  D02fKR7>y 

1  (D%k)2RT/i  RT 

With  Eq.  (20),  the  pressure  distribution  in  the  porous  cathode  can 
be  analytically  calculated  without  any  numerical  techniques. 

Fig.  4(a)  shows  the  distribution  of  pressure  in  the  SOFC  cathode, 
calculated  by  numerical  method  (DGM)  and  theoretical  solutions 
(Eq.  (20)).  It  is  found  that  the  pressure  distribution  calculated  by 
the  theoretical  method  (Eq.  (20))  agrees  very  well  with  that  by 
numerical  method  (DGM).  Thus,  the  theoretical  solution  devel¬ 
oped  in  this  paper  is  useful  to  calculate  the  pressure  distribution 
quickly  and  accurately.  As  expected,  the  pressure  decreases  along 
the  cathode  depth  and  the  pressure  gradient  increases  with  increas¬ 
ing  current  density  due  to  electrochemical  reduction  of  02  at  the 
TPB.  Obviously,  this  pressure  gradient  is  beneficial  for  02  perme¬ 
ation,  leading  to  higher  02  partial  pressure  in  the  porous  cathode 
(Fig.  4(b))  and  lower  concentration  overpotentials  (Fig.  5(a)).  The 
difference  between  the  two  cases  is  also  considerable.  For  example, 
at  a  typical  current  density  of  5000  Am-2,  the  cathode  concentra¬ 
tion  overpotential  considering  pressure  gradient  is  about  0.0049  V, 
while  it  is  0.0054  V  (10%  difference)  for  the  case  with  the  pressure 
gradient  neglected  (Fig.  5(a)).  Different  from  the  anode,  the  cathode 
concentration  overpotential  increases  with  increasing  temperature 
(Fig.  5(b)).  This  is  because  of  lower  gas  density  at  higher  tempera¬ 
ture  (Eq.  (2)),  despite  the  fact  that  the  effective  diffusion  coefficients 
increase  with  increasing  temperature  (Eqs.  (5)  and  (6)).  In  addition, 
the  difference  between  the  two  cases  is  found  independent  of  the 
operating  temperature. 

4.  Conclusions 

A  comparative  study  has  been  conducted  to  study  the  impor¬ 
tance  of  pressure  gradient  in  modeling  of  SOFC  concentration 
overpotentials.  It  is  found  that  the  pressure  gradient  in  the  porous 
anode  of  an  SOFC  running  on  CH4  is  significant.  This  pressure  gradi¬ 
ent  impedes  the  transport  of  H2  to  the  TPB  but  facilitates  transport 
of  H20  from  the  TPB  to  the  outer  layer  of  the  anode.  As  a  result, 
neglecting  this  pressure  gradient  can  underestimate  the  concentra¬ 
tion  overpotential  by  about  20%  under  typical  working  conditions. 
This  underestimation  would  become  more  serious  at  a  higher  tem¬ 
perature.  It  is  thus  concluded  that  the  pressure  gradient  in  the 
porous  anode  of  an  SOFC  running  on  hydrocarbon  fuels  with  inter¬ 
nal  reforming  must  be  considered.  At  the  cathode,  neglecting  the 
pressure  gradient  can  overestimate  the  concentration  overpotential 
by  about  10%  under  typical  working  conditions. 
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